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Abstract 

We demonstrate that the boundary of two optical isomers with opposite 
directions of the gyration vectors (both parallel to boundary) can support 
propagation of electromagnetic wave in the direction perpendicular to the gy- 
ration axes (Cotton-Mouton geometry). The components of electromagnetic 
field in this wave decay exponentially into both media. The characteristic de- 
cay length is of the order of the Faraday rotation length for the propagation 
along the gyration axis. The remarkable property of the boundary wave is its 
chirality. Namely, the wave can propagate only in one direction determined 
by the relative sign of non-diagonal components of the dielectric tensor in con- 
tacting media. We find the dispersion law of the boundary wave for the cases 
of abrupt and smooth boundaries. We also study the effect of asymmetry be- 
tween the contacting media on the boundary wave and generalize the results 
to the case of two parallel boundaries. Finally we consider the arrangement 
when the boundaries form a random network. We argue that at a point, when 
this network percolates, the corresponding boundary waves undergo quantum 

derealization transition, similar to the quantum Hall transition. 
PACS numbers: 72.15.Rn, 73.40.Hm, 78.20.Ek 
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Optical properties of gyrotropic and optically active media are well known and described 
in many textbooks (see e.g. Ref. |1|). A homogeneous gyrotropic medium is characterized by 
the following relation between the displacement vector and electric field D = eE, where the 
tensor e has the following form 

e ig 

-ig e • (1) 
\ e lj 

For a wave propagating along the z-direction, the wave equation 

V(VE)-V 2 E = -i^ (2) 

has two circular polarized solutions characterized by the dispersion laws 

to 



k± = -V £ o ± 9- (3) 
c 

A linear polarized wave, incident on this medium along the axis of gyration, z, splits into 
two circular polarized waves propagating with velocities cj (so±g) l l 2 . Since g is always much 
smaller than Eq, this can be viewed as rotation of the plane of polarization (Faraday effect). 
The distance at which the plane of polarization is rotated by 90° is equal to l u = iiy/e c/guj, 
where uo is the frequency of the wave. The direction of rotation is determined by the sign 
of the nondiagonal component, g, of the tensor ([!]). Correspondingly, in the case of optical 
activity, there exist left-rotating (levorotatory) and right-rotating (dextrorotatory) media. 
Two modifications of a crystal differing by the sign of g are called optical isomers. For 
example, the parameters of tensor e for the quartz crystal are known to b J e = 2.3839, 
£\ = 2.4118, g = 1.1 ■ 10 -4 . Consequently, l u ~ 4mm for yellow light (A = 589nm). Isomers 
with tensor e in the form (1) can only belong to the certain crystalline symmetry classes 
listed e.g. in Ref. ||. 

When the direction of propagation is perpendicular to the gyration axis, there also exist 
two types of solutions of the wave equation (fj). Firstly, there is a trivial solution for which 
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only E z component is nonzero and the spectrum is u = kc/ 'y/e 1 . The nontrivial solution 
(Cotton-Mouton effect) corresponds to polarization perpendicular to the axis of gyration 

E=*°( vtl7C> )e " (4) 

with the dispersion law 

kl + k 2 y = ^(e -f). (5) 

c to 

The phenomenon of surface electromagnetic waves is also a well-studied subject (see 
e.g. Ref. 0). Surface waves can propagate along the boundary of two isotropic media with 
dielectric constants e a {uo) and £&(u;) when two conditions are meti: (i) e a {uo) ■ Sb{<jo) < and 
(ii) e a (uj) + Et,(uj) < 0. The polarization of these waves is TM (transverse-magnetic). For 
such waves the localization of electromagnetic field at the boundary is provided by negative 
sign of e in one of the contacting media. Interesting examples of surface electromagnetic 
wave at the boundary of two transparent media were given in Refs. |] and |5|. In Ref. f| the 
boundary between isotropic and uniaxial media is considered. In Ref. [| both contacting 
media represent uniaxial crystals with certain mismatch in the directions of optical axes. 

In the present paper we show that the boundary of two optical isomers with opposite 
directions of axes of optical activity (both parallel to the boundary) can also support propa- 
gation of localized electromagnetic waves. In contrast to the conventional surface plasmonsi, 
these states exist for positive values of components e$ and E\. The distinguishing feature of 
these boundary waves is a unidirectional character of propagation, i.e. they propagate along 
the boundary of two media with opposite gyration vectors only in one direction which is 
determined by the relative sign of g. 

Let the plane x = be a boundary between two isomers. For generality we consider g 
to be some function of x , g(x), changing its value from g(x) = go to g(x) = —g within 
some region around x = 0. In the Cotton-Mouton geometry, when only x and y components 
of electric field are nonzero, the system of equations for E x , E y which follows from Eq. (^) 
takes the form 



d 2 E,, d 2 E T 



OJ 



dxdy dy 2 
d 2 E x d 2 E,„ 



OJ 



(e E x + ig{x)E y 



{e E y - ig(x)E x ). 



(6) 
(7) 



dxdy dx 2 c 

Let us search for a solution propagating along the y axis: E x oc e lky , E y oc e lky . Then from 
Eq. (6) we can express E x through E y 

iui 2 E y g{x) 



E r 



ikc 2 ^ 

ox 



e oj 2 — k 2 c 2 



(8) 



Substituting Eq. (|8|) into Eq. (7) we get 



d 2 E y 
dx 2 



+ 



k dg(x) 
e dx 



+ 



c 2 e 



(9 2 (x)-g 2 ) 



E„ 



C 2 Eq 



E„. 



(9) 



Note now, that Eq. @ has the form of the Schroedinger equation with an effective potential 



U{x) 



k dg(x) 

Eq Ox 



+ 



C 2 E 



(10) 



and energy 
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k 2 . 



(11) 



Consider first the simplest case when the transitional region between two media is in- 
finitely narrow: g(x) = —g (29(x) — 1), where g > and 9(x) is a step-function (see Fig. |TJ). 
Then U(x) takes the form U(x) = —2g k5(x)/EQ. It is seen that the potential is attractive 
for the positive values of k and repulsive for the negative k. This illustrates the above state- 
ment about the chirality of the boundary states. It is well known that there exists only a 
single bound state in attractive 5-potential with an arbitrary magnitudei. The correspond- 
ing "binding energy" is equal to S = —g 2 k 2 jE\. This leads to the following dispersion law 
for electromagnetic wave propagating along the abrupt boundary between two media with 

g = ±go 

ck 



OJ 



(12) 



The origin of the localized boundary state Eq. (|TJ) is the following. It is seen from Eq. ([S]) 
that the maximal possible wave vector inside each of the contacting media is equal to 



*c ( £ o~ f^) 1 ^ 2, ^ e wave vector determined by Eq. flI2|) exceeds this maximal wave vector. 
As a result, both components of electric field decay to the left and to the right from the 
boundary x = 

E y (x) = E y (0)e~^ (13) 

' t 2 2 \ 

E x {x) = Sl gn(x p £ ° + go) E y (0)e-^, 

where q = \kgo\/e (Fig. [I]). The characteristic decay length, 1/q, can be expressed through 
the length of rotation of plane of polarization of light in the Faraday geometry: 1/q = l u /iv. 
It is easy to establish certain properties of chiral boundary states. 

a) Asymmetry between the contacting media. Suppose that the mirror symmetry between 
the contacting media is lifted, say, due to external magnetic field applied in the z direction. 
This amounts to the following modification of g(x) : g(x) = —go{29{x) — l)+gi, where g\ is 
proportional to the magnetic field. Then it is straightforward to check that the dispersion law 
Eq. ([Tj|) remains unchanged, while the magnitude of g± remains smaller than go- When g\ 
exceeds g the bound state disappears abruptly. Despite the dispersion law being insensitive 
to the asymmetry for g x < g , the distribution E y (x) changes strongly. For g x ^ instead 
of Eq. ( Jf3| ) we have 

E y (x) = E y (0)exp\ — (go-gi)x\, x>0 (14) 



E y (x) = Ey(0)exp\——(g + gx)x], x<0 (15) 



b ) Two boundaries. Consider now the case of two boundaries 

g{x) = -g , \x\ > 1/2 (16) 
g{x) = g , \x\ < 1/2 (17) 

It is straightforward to analyze Eq. @ in this case. Indeed, the effective potential takes the 
form 



U{x) = 2g k[5(x - 1/2) - 6(x + 1/2)}. 
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(18) 



The corresponding dispersion law for the boundary state becomes 

C 2 E 



(19) 



where q\ is determined by the equation 

1 



Qi 



9l 



-2 gi l 



^9o 
c 2 e 



1 _ g l 

F 2 



(20) 



In the limit I — > oo we return to the dispersion law Eq. (J12f) , which for positive k corresponds 
to the wave propagating along the boundary x = 1/2; negative k corresponds to the wave 
propagating along the boundary x = —1/2. Inspection of Eq. (|20|) shows that it has solution 
for an arbitrary small I. In the case when the distance between boundaries is much smaller 
then the localization length we get from Eq. (f20|) 

/ 



qi = ^ <<l: 



(21) 



This means that both states (with positive and negative k) are "weakly bound". This is 
illustrated in Fig. ^| for different l/l^ and positive k. For negative k the distribution of 
electric field corresponds to the change x — > —x. 

c) Two boundaries with external magnetic field. We have also studied the suppression of 
the states associated with a pair of boundaries by an external magnetic field. In this case 
the critical value of g\ depends on the distance between boundaries. The electric field in the 
region \x\ < I represents the superposition of exponents exp( 2 y^) and exp(— z y^), where the 
dimensionless parameter k satisfies the equation 



K 



i | 2 9l 2g 2 Q l + exp(-2LV) / | Ag x 



k 2 g 



Eq 1 — exp(— 2L*k) 





= 2 


V 9o* J 





(go - gi) 



(22) 



with L* = rrl/l^. In Fig. ^| we present the critical line of gi/go as a function of L*, which 
determines the domain where solution of Eq. (^) exists. Note, that the boundaries are 
asymmetric with respect to the change of the direction of magnetic field (the sign of #i). 

d) Smooth boundary. Next we consider the case when the boundary is smooth and has a 
characteristic width of b. This, for example, can be a result of mutual diffusion of isomers. 
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Then g(x) can be modeled by g(x) = — go tanh(x/6). Thus, for the effective potential U(x) 
we get 



U{x) 



g (kc 2 + u 2 g b) 



e bc 2 cosh 2 (x/b) 

The solutions of the Schroedinger equation with potential Eq. (|2~3"D are well known@ 



(23) 



£n ~ 46 2 



l + 4(^ + ^) -(2n + l) 



1 2 



(24) 



e e c* 

The mode with n = has no threshold frequency. The threshold frequencies, u n , for the 
modes with n > are determined from the condition £„ = 



1 



(2n+l) 2 -4 



\ 



1 _ ^ 

-0 



(25) 



The corresponding dispersion laws, uj n (k), can be conveniently presented after introducing 
a dimensionless frequency, Q, and the wave vector Q 



U = — uj, 



Q = rbk, 



(26) 



where the dimensionless parameter r is defined as 



T 



go 

so' 



(27) 



Then from Eqs. ( |TTD and ( 21 ) we have 

n(Q) = I 



(2n + 1)t 2 + J(2Q + 1 - r 2 ) 2 - An(n + l)r 2 (l - r 2 ) 



1/2 



-(1 + 4Q) 



(28) 



The dispersion law for the first five modes is shown in Fig. |j. 

Let us consider qualitatively the situation when nondiagonal component of the tensor e is 
a random function of both coordinates x and y. Suppose for simplicity, that the correlation 
length of g(x, y) is much bigger than Z w . Then it is obvious from the above consideration 
that the boundary waves would circulate along the contours g(x,y) = 0. If the average 
value g{x,y) is negative (Fig. 5a) or positive (Fig. 5b), and comaparable to (g 2 {x,y)J ^ , 
these contours are disconnected, and, correspondingly, the boundary waves are localized. 



As g(x,y) approaches zero, the contours, defined by the condition g(x,y) = 0, grow in size 
and form a network. The correlation size of the network divergies as (g(x,y)^J \ where 
v\ pa 4/3 is the critical exponent of the percolation problem in two dimensions^. Along 



with decreasing of g(x, y), the neighboring contours g(x, y) = come closer than l w , and the 
coupling of the boundary waves, encircling these contours, becomes increasingly important 
(above we studied this coupling for the simplest case of parallel boundaries). As a result 
of this coupling, the interference of different unidirectional paths comes into play. The 
crucial role of interference effects, allowed by the coupling, was first pointed out in a seminal 
paper Ref. || in relation to the integer quantum Hall effect. In Ref. || unidirectional waves 
modeled the motion of a two-dimentional electron in a strong perpendicular magnetic field 
and a smooth random potential (edge states). It was demonstrated in Ref. § that, with 
interference taken into account, the derealization transition in the system of edge states 
occurs at discrete energies, i.e. it remains infinitely sharp (as it is the case for the classical 
percolation). However, due to the coupling and interference, the size of the eigenstates in 
the critical region (localization radius) is much bigger than the correlation radius of the 
network and divergies with the exponent^ v-i ~ 2.3. The correspondence between the edge 
states of electrons and the boundary electromagnetic waves allows us to conclude that with 

-v 2 



g(x,y) — > the radius of localized boundary waves behaves as (g(x,y) 

By analogy to the integer quantum Hall transition, which originates from the competition 
of unidirectional motion of an electron in a magnetic field and quantum interferenceiE, the 



derealization transition of the boundary waves at g(x, y) = (when optical activity of the 
system is zero on average) can be called quantum Cotton-Mouton effect. 

Note in conclusion, that throughout the paper we assumed that the system is uniform in 
the z-direction. The situation relevant for experiment is a thin-film geometry, in which two 
contacting optically active media are confined within the region \z\ < | with the thickness, 
d, of the film being much smaller than l u . If Eq exceeds the dielectric constants of the 
media, between which the film is sandwiched, the solutions of the Maxwell equations are 
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the waveguide modes, propagating along the plane of the film. The components of electric 
and magnetic fields in these modes are confined within the region of the order of d. Suppose 
first, that g is constant within the film: g(z) = <?o$(f — \z\). To modify the Cotton-Mouton 
dispersion law Eq. (|5|) to the case of the waveguide mode propagating in the y-direction, it 
is convenient to rewrite the system of Maxwell's equations in the following form 



d 2 E x 
dz 2 



(e (z)— - k 2 )E x 



uj 2 g 2 (z) u g(z) dB x 



-E 



c 2 Eq{z) x c s (z) dz ' 



(29) 



d ( 1 8B X 



e {zf--k 2£ MV l 



uj d 



g{z)E x 



(30) 



dz \sq(z) dz J \ c 2 J c dz 

where Eq(z) and £\(z) describe the profile of the diagonal components of e in z-direction. 
For g = the right-hand sides in Eqs. ([2^,[3(]) are zeros, so that the above equations yield 
the sets of TE and TM waveguide modes, respectively. With the right-hand sides included, 
the correction to the wave vector k is quadratic in g. For TE mode with a number n after 
some algebra one can get the following dispersion law 



k 2 = {k^)) 2 -^[{g^f 



(31) 



where 



and 



si 



rd/2 , (4 n) )2 

J-d/2 az e (z) 



Poodz(E^) 2 



I-oo dz 



R (m) „(n) 

eo(z) 



? IZ,dzB x m) E x n %(e (z)-y 2 ) 



f d / 2 J~ E™ &Bi m) 
-d/2 uz £o (z) dz 



(32) 



rd/2 , F(n) d (b^\ 

J_ d/2 dzti x ' dz { £o(z) ) 



IZodziE^J^dz^^ 



£o(z) 



(33) 



(2) N2 



The correction (g e L) 2 originates from the g 2 -term in the r.h.s. of Eq. (p9[) , whereas (g e ff 
results from the mixing of TE mode n with all TM modes. Generally speaking, g^L and 
g e jj are of the same order. This means that, when g depends on x and changes sign, 
the corresponding boundary wave would represent a mixture of TE and TM modes. The 



situation simplifies if, for numerical reasons, g^L appears to be much bigger than g e Jt- In 
this case the dispersion law for the boundary wave, analogous to Eq. flT2]), takes the simple 
form k = kj^(u), and the electromagnetic field decays away from the boundary x = as 
exp(-^ ) / |x|). 
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FIGURES 

FIG. 1. The boundary of two optical isomers (x = 0). Boundary wave propagates along the 
y-axis. The decay of E x (dashed curve) and E y (solid curve) components of electric field away 
from the boundary x = is shown schematically. 

FIG. 2. E y {x) component of electric field for the case of two boundaries. The location of the 
boundaries are x/l w = ±1.25 (dotted line); x/l^ = ±0.75 (dashed line); x/l w = ±0.3 (solid line). 

FIG. 3. Critical line gi/go as a function of L* for two opposite directions of external magnetic 
field. 

FIG. 4. Dispersion relations for the guided modes with n = 0,1,2,3,4 are plotted using 
Eq. (28). For illustration purposes we chose unrealistically large value of parameter r = 0.5. 

FIG. 5. The contours g(x, y) = are shown schematically for the case: (a) g(x, y) < and (b) 
g(x, y) > 0. The arrows show the direction of propagation of boundary waves. 
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